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Abstract 

We study the production of heavy quarks in deep-inelastic scattering within perturbative QCD. 
As a novelty, we employ for the first time the running mass definition in the MS scheme for 
deep-inelastic charm and bottom production. We observe an improved stability of the perturba- 
tive expansion and a reduced theoretical uncertainty due to variations of the renormalization and 
factorization scales. As our best estimate we extract from a global fit to fixed-target and HERA 
collider data for the charm-quark an MS mass of m c {m c ) = 1.01 + 0.09 (exp) + 0.03 (th) GeV. 
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1 Introduction 



The production of heavy quarks in deep-inelastic scattering (DIS) is an important reaction and has 
been measured with high accuracy in several fixed-target experiments and at the HERA collider. 
Within perturbative QCD, the production of charm and bottom quarks proceeds in neutral (NC) 
or charged current (CC) reactions via lepton-parton scattering and the exchange of a virtual boson 
y*/Z or W ± with space-like momentum. A detailed understanding of the production mechanism 
sheds light on the underlying parton dynamics in QCD. In the LHC era this, perhaps, is the most 
important aspect, because DIS heavy-quark production provides core constraints in global fits on 
the parton distribution functions (PDFs) even at the terascale. 

Thus, it is of paramount importance to provide precision predictions which, of course, have to 
rely on higher order radiative corrections. Our theory predictions for heavy-quark production in- 
clude next- to-leading order (NLO) QCD corrections DTH21 and, in the case of NC DIS, even partial 
information at next-to-next-to-leading order (NNLO), which comprises in particular all logarith- 
mically enhanced terms near threshold 01, and all explicit dependence on the renormalization and 
factorization scales. 

However, precision predictions, must also address the uncertainty due to the non-perturbative 
parameters, such as the aforementioned PDFs, the value of the strong coupling constant a s and 
the mass m of the heavy quarks charm, bottom and top. It is precisely with respect to the latter 
aspect, that we wish to improve the current state-of-the-art. Namely, we employ the short-distance 
(so-called MS) mass in our treatment of heavy-quark DIS. In this manner we provide a crucial 
link, which has long been missing, for the comparison of the heavy-quark masses entering in DIS 
and the determination of PDFs in global fits on the one hand, and, on the other, those obtained 
from other determinations, e.g. in e + e~-collisions or by means of lattice computations. 

Traditionally, perturbative corrections to hard scattering processes at hadron colliders have used 
the so-called pole mass of the heavy quark as a definition of the mass parameter. The pole mass 
is popular, because it is well defined at each finite order of perturbation theory and it is introduced 
in a gauge invariant way. However, as is well-known since long, the concept of the pole mass in 
QCD has intrinsic theoretical limitations. Because of confinement no free colored quarks exist, i.e. 
they do not appear as an asymptotic state of the S -matrix. It has been shown that the use of the 
pole mass leads to a poorly behaved perturbative series 0, because observables in hard scattering 
processes become sensitive to momentum regions of the order of the QCD scale Aqcd- 

Alternative mass definitions offer a solution to this problem. The most prominent example is 
the MS mass m(// r ), which is to be evaluated at the (renormalization) scale fi r , where n r » Aqcd, 
and which is free of ambiguities of order Aqcd- For inclusive cross sections at short distances the 
appropriate scale choice for the running mass m(ju r ) is /u r = m, where the renormalization group 
evolution for the scale dependence of the mass converges even for scales as low as the charm-quark 
mass. As a benefit of theory predictions using the MS mass one observes an improved stability of 
the perturbative series with respect to scale variations as compared to the result in the pole mass 
scheme. 

In this study, we start off with QCD predictions for heavy-quark DIS at NLO Q3-[3|| and approx- 
imate NNLO [4], which have been computed with a pole mass. Subsequently, we improve the per- 
turbative description by converting from the pole mass scheme to the MS scheme (see Refs. (6]-[8]| 
and references therein). The necessary scheme transformation follows closely similar recent work 
for top-quark production at hadron colliders, cf. [OiTOj for implementation details. From a global 
fit of the parton distribution functions to fixed-target and HERA collider data we extract for the first 
time an MS mass for the charm-quark m c (m c ). Our best estimate for m c (m c ) is consistent with the 
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world average [fTTTl within the quoted range of errors. Moreover, PDFs determinations can benefit 
from reduced uncertainties due to precise heavy-quark masses for charm and bottom, m c and nib, 
and we comment on the implications for W ± and Z gauge boson production at the LHC. 



2 Heavy-quark DIS in perturbative QCD 

For NC DIS pair-production at leading order (LO) proceeds through photon-gluon fusion as 

g(p)+7*(q) -^qi+qi, (l) 

which is a 2 — » 2 process starting off at order a s in QCD and, of course, involving the overall power 
a for the QED coupling and the quark fractional charges. 

For CC DIS, on the contrary, heavy-quark production at parton level proceeds in Born approx- 
imation in a 2 — » 1 reaction as 

qi(p) + W*(q)^q 2 , (2) 

where the initial quark q\ is light, the final state quark q 2 is heavy and the coupling to the W-boson 
involves the usual parameters of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. 

The well-known kinematical variables are Bjorken x and Q 2 defined by the momenta p and q 
of the incoming parton and the off-shell boson, 

O 2 

Q 2 = -q 2 >0 X = T—^ (3) 

2p-q 

and the cross sections is conveniently parametrized in terms of the heavy-quark DIS structure 
functions Fk, k = 1,2,3, which depend on x, Q 2 and the heavy-quark mass m. In the standard fac- 
torization approach to perturbative QCD the structure functions Fk can be written as a convolution 
of PDFs and coefficient functions, 

..max 

F k (x,Q 2 ,m 2 )= £ f jfi^fjCk^Z^lfij), (4) 

i=q,M x 

where the renormalization scale is denoted as p r and the PDFs for the parton of flavor i at the 
factorization scale p , as fj(x,p 2 ). Depending on the kinematics in Eqs. (QQ), (O, the integration 

range over the parton momentum fraction z extends to z™q = 1/(1 +m 2 /Q 2 ) or = 1/(1 + 
Am 2 IQ 2 ). The kinematical variable £ in Eq. © is given as 

Q 2 

m z 

The coefficient functions Ck,i of the hard parton scattering process in Eq. © can be computed 
in a perturbative expansion in the strong coupling constant a s = a s (ju r ). Currently, we have for 
both cases, NC and CC, the complete NLO corrections available with full dependence on the 
heavy-quark mass m, see Refs. [QQ| and OEl, which we use in our description of the heavy-quark 
structure functions Fk- Specifically, in the NC case, we use the code of Ref. [fT2l (see [fT3l for 
minor corrections). 

Beyond NLO, partial results are available, although the complete NNLO corrections are not 
known to date. In the asymptotic limit m 2 /Q 2 —* fully analytic results have obtained at NLO, 
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see [ l4] - fT6l and at NNLO for the lowest even-integer Mellin moments IfTTl . For parton energies 
close to the production threshold, s - Am 2 , soft gluon improvements at NNLO are long known to 
be important [fT8l[T9l . For NC heavy-quark production the convolution of the coefficient functions 
(especially C2,g) and the gluon density is dominated by rather low partonic of-mass energies s and 
the corresponding soft logarithms in /? = (1 -4m 2 /s) 1 ^ 2 at NNLO have recently been completely 
determined flU. In our description of the NC structure functions F k (which supersedes our earlier 
studies |fT9l ) we include these latest improvements [0 together with complete dependence on the 
renormalization and factorization scales, see e.g. BUdSl. This approximation to NC DIS we call 
NNLO a pp r ox- Soft logarithms have also been studied for the CC case, see GUI . However, in the 
kinematical range of the currently available CC DIS data, they are numerically less important and 
we do not include them here. Hence, for CC DIS, we confine ourselves to NLO accuracy only. 

The mass parameter in the structure functions F k in Eq. © is the pole mass of the heavy 
quark, which requires m to coincide with the pole of the heavy-quark propagator at each finite 
order in perturbation theory. In this way, that value of the mass itself is strongly depended on 
the perturbative order. Moreover, it has intrinsic uncertainties of order Aqcd/'^- The perturbative 
description of heavy-quark DIS can be improved, however, by performing a scheme change from 
the pole mass to the MS scheme. 

The starting point of this conversion is the well-known relation between the pole mass m and 
the running mass m(ju r ) in the MS scheme 

m = m(n r ){\ +a s (ju r )d (r) (M r ) + a s (ju r ) 2 d (2) (ju r ) + ...), (6) 

where the coefficients <i® of the perturbative expansion in a s are actually known to three-loop 
order HHH. 

Let us start with the NC case. We will derive explicit formulae through NNLO for the depen- 
dence of the structure functions on the MS mass m(m). In doing so, we follow similar recent work 
for the pair-production of top-quarks at hadron colliders [0[TOj. For the pole mass m we have 
(suppressing all other arguments), 

F k (m) = a s Ff\m) + a 2 F%\m) + a s 3 Ff\m) , (7) 

which we can convert with Eq. © to the MS mass m(m) (for simplicity abbreviated as m) according 
to 

F k (m) = a s Ff\m) (8) 



+a 2 (if } (m) + md m d m Ff\m) ) 

\ " m=m) 

+a s 3 ^Ff\m) + md {2) d m Ff\m) 



+ md w d m F ( ^\m) 



+\{md^) 2 d 2 m Ff\m) 



where the coefficients d^> have to be evaluated for /u r = m (corresponding to the scale of a s ). 

In the NC case, the coefficient functions in Eq. © have a perturbative expansion to NNLO in 
the strong coupling a s = a s (pi r ), 

C k MzU,»l» 2 f) = a s cf) + a 2 C { ^ + a s 3 eg , (9) 
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where 77 denotes the distance to partonic threshold. The partonic center-of-mass energy reads 
s = Q 2 (l/z-l), so that 

s 

Am 2 Am 2 \ : 

Mass dependence resides in the coefficient functions implicitly in 77 and £ as well as in the 

>}>' 



mz) = -^-i = -^(--i\-i. do) 



factorization scale dependent part (commonly appearing through the ratio /j. 2 /m 2 ). Thus 



d m Ckj(r](z),^,id 2 ,fi 2 f ) = (d m rj) d^Ckj + (d m £) dfC k j - —d^Ckj ■ (1 1) 

In this way the explicit expression for the first-order derivative d m in Eq. © becomes 

wfw - z /f^fK^^-^-^s} <12) 

i=q,q,g x 

+ (d m z™) X -fil-)^ , 

Z \Zl ' z =z mm 

where we have suppressed all arguments in the coefficient functions Ck,i for brevity. The derivative 
d n has been turned into the derivative d z thanks to Eq. (flOl) and all partial derivatives in Eq. (ITT]) 
have been made explicit. 

The boundary term in Eq. (Tl2l) vanishes explicitly for the NLO scheme transformation in 
Eq. ©, i.e. for d m F^\m), because the Born contribution behaves for fi — » as C^ g ~ 0(J3) and 

cj.°> ~ 0(J3 3 ). Thus, cf ] g vanishes in the last line of Eq. (TJH) if evaluated at z = z max We note 
however, that the boundary term may be completely removed to all orders by means of partial 
integration with respect to the PDFs. With integration-by-parts in z we find for Eq. (fT2l) . 

w» - 1 Jf>i(f){K-5fvs-^} 



^max 

2 JtM*-**) 



* X 

The NNLO scheme transformation in Eq. ® requires as the only additional ingredient the 

-(0- 

k 



second derivative d 2 for the Born term f{ (m), i.e. 



£ff(») = -4 ?VfW + -^& w Ff(m)L (14) 



/;/- 



1 V 2 K J m 



which can be quickly evaluated using Eq. (TOT) and computing the explicit derivative of the coef- 
ficient functions cfK We note, that d^cf g ~ Oifi) and d^cf g ~ OQ3 3 ), so that they also vanish if 
evaluated at z = z max - 

Thus, we are finally in a position to put everything together through NNLO and we arrive at 
the following explicit expression for Eq. ©, 

F k = a,Ff + a*F?>+a*Ff> (15) 
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+- 2 .E / 7^(f)^{=«2 



'=q,q,g 



+a x 



+a x 



+a x 



X 

„max 



2 JiM'-n^H*®*-* 



i=q,q,g 



2d m C {l) k 

i,k 



x 

jnax 



E / f [<)M*}+fjr(f)M 2 Wi^}- 



•=q,q,g 



where the renormalization scale has been fixed at = m(m), i.e. a s = a s (m(m)) for the strong 
coupling constant. The full dependence on \x r can be constructed using the renormalization group 
equation. With the standard expression for the running coupling (and the coefficients of the QCD 
beta function), it is easy to restore the complete renormalization scale dependence of a s in Eq. (fT3T) . 

a s (m(m)) = a s (p. r )(l + a s (^ r )L^ + a s Qu r ) 2 (J3iL R +/3lL 2 R )^ , 



(16) 



where we have abbreviated L R = \n(ji 2 /m(m) 2 ). 

As explained above, Eq. (fT5l) is exact to NLO. At the NNLO level, cl : is currently unknown. 
Our approximation NNLO ap prox in Eq. (TT5l) uses the threshold enhanced terms of Ref. flU to esti- 
mate the dominant corrections at two loops. All explicit dependence on the renormalization and 
factorization scales and the terms accounting for the scheme transformation from the pole to the 
running mass at NNLO are, however, exact in Eq. (fTSb . 

The CC case is conceptually much simpler. Moreover, in our analysis, we confine ourselves to 
the NLO case only so that the relevant formulae is much shorter. One commonly defines structure 
functions Tk which are related to the FkS via the following relations: 



F\ = Tx =a l Ti, 
F 2 = 2xT2 = a 2 T2, 
F 3 = 2T3 = a 3 T3, 

i.e. a\ = 1, «2 = 2^ and a 3 = 2 and where we have introduced the quantity 



(17) 
(18) 
(19) 



X = T 

so that Bjorken x varies in the range < x < A and 

1 



A = 



1 + wL 1 +£ 
1 + e 2 



(20) 



(21) 



In the CC case, the coefficient functions corresponding to Eq. © can be expanded to NLO in 
the strong coupling with a s = a s (jji r ), 



C k>i (z,^r^f) = CfJ + a s C^.. 



(22) 
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where cf } q ^S(l-z) (up to the CKM parameters) and cf ] g = for k = 1,2,3 due to Eq. ©. The 

expressions for cP are all given in Refs. (2][3]]. As in the NC case, the conversion to the MS mass 

starts off from a perturbative expansion as in Eq. ©, which we can convert to the MS mass as in 
Eq. ©. Restricting ourselves to NLO, we relate 

F k (m) = Ff ) (m) + a s F ( k l \m), (23) 

to 

F k (m) = Ff\m) + a s [F^\m) + md {l) d m Ff\m) V (24) 

\ m=m) 

where matching is done again at the scale fi r = m and we just have to provide the explicit expression 



md (l) d m akF\ \m) 



= 2(x-x)d (1) J] \Vi q2 \ 2 {d x a k fi(x)) 



(25) 



where the sum ranges over all light flavors q\ , while q2 is heavy. V qi qi denotes the respective CKM 
matrix element and the coefficients a k are defined in Eqs. (fT71)-(fT9l). The partial derivative from 
dmX has again been made explicit. 

This completes our discussion on the explicit conversion of the heavy-quark DIS structure 
functions from the pole to the MS mass. Note, that we have confined ourselves entirely to the 
so-called fixed- flavor-number scheme (FENS), i.e. to a situation, where we work with a fixed 
number nt of light-quark flavors. This is an absolutely adequate approach e.g. for the analysis of 
existing DIS data. In contrast, there exist variable flavor number schemes (VFNS) which relate 
the DIS structure functions for n t light flavors to those for rif + 1 light flavors (see [|2Tll22l for an 
extensive discussion). In a VFNS the necessary matching involves certain massive operator matrix 
elements which, for consistency also need to be evaluated in the renormalization scheme with a 
running mass. To NNLO, all relevant formulae has been given in Ref IfPTTl . Also, in practice, one 
often relates the matching scale of the strong coupling to the heavy-quark mass, i.e. a 
a/- n f \ian) with some constant k. Then the necessary decoupling coefficients depend beyond 
NLO on the chosen mass renormalization (see e.g. |j23l). In summary, the implementation of a 
running mass for heavy-quark DIS in a VFNS can be performed in a straight forward manner and 
poses no further problems. 



3 Results 

We are now in a position to look at the phenomenological implications of the MS mass and also 
to discuss the role of the heavy-quark mass parameter in PDF fits. Despite of its short-comings, 
all current global fits of PDFs employ the pole mass scheme for the heavy quarks. To illustrate 
this point, we summarize in Tab. \T\ the values taken by the six groups which are currently ac- 
tive in global fits PDF: ABKM [EH, HERAPDF (21, GJR (251, MSTW (261, CTEQ [27] and 
NNPDF [28]. All numerical values, especially those for the charm quark mass, are systematically 
lower than the pole masses obtained from the particle data group (PDG) values for the world aver- 
age [1 1J. Simple kinematical considerations show, that a smaller charm quark mass can potentially 
compensate large missing higher order perturbative corrections. The latter have been shown to be 
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[GeV] 


PDG 


ABKM 


GJR 


HERAPDF 


MSTW 


CTEQ 


NNPDF 


m c 
m b 


1 66 +009 

luu -0.15 

4 7Q +0.19 
*-' y -0.08 


i <r +0.25 
J -0.25 

4 5 +0-5 
-0.5 


1.3 
4.2 


1.4 
4.75 


1.4 
4.75 


1.3 
4.5 


V2 
4.3 



Table 1: The pole mass values taken as input in recent global fits of PDFs. The quoted PDG values are 
obtained from the MS values in Eqs. (f26T >. (1271 ) using the two-loop conversion in Eq. ©. 

sizable at NLO and even at NNLO ap prox when using a pole mass. Evidently, global fits of PDFs 
which incorporate heavy-quark DIS data are very sensitive to the theory treatment of heavy quarks 
and depending on the chosen mass parameter (and scheme, see ETTl for a discussion of FFNS and 
VFNS) the resulting differences in the PDFs can be sizable. 

In Figs. [Hand [2] we study first the case of charm quark electro-production, i.e. electron-proton 
scattering. We plot the NC charm structure function F p at LO, NLO and NNLO ap prox using the 
3-flavor PDF set of Ref. ETTl . We have been careful to restrict the kinematics in x and Q 2 (in 
Figs.IHEle 2 = lOGeV 2 , x = 1(T 3 ) to the region, where our threshold approximation underlying 
the NNLOapprox prediction is under control, see H18[|19H . 

Comparing the central values of the predictions for F% as a function of the pole mass in Fig. Q] 

(left) with those for an MS mass m c (m c ) in Fig. \T\ (right), we observe a much improved convergence 
of the perturbative expansion in the latter case. Already at NLO the size of the QCD corrections 
is much reduced for an MS mass. E.g. for a pole mass of m c = 1.5 GeV we find for F p a relative 
increase of 32% at NLO over the LO prediction and another 13% for the NNLO ap prox prediction 
normalized to the NLO one. This is to be compared with the numbers for the MS mass. At m c (m c ) = 
1.3 GeV we find relative corrections of 17% at NLO and only 6% at NNLO ap prox- Likewise, as 
we vary the renormalization and the factorization scale fj. r and fj.^ independently by a factor of two 

around a central value chosen to be ji 2 = ji 2 = Q 2 + 4m 2 in Figs. [Hand [2] we note a substantial 
reduction in the spread to the predictions and a greatly reduced theoretical uncertainty when using 
the running mass. This scale variation is illustrated by the respective bands in Fig.[2l 

Similar observations hold also for the NC DIS production of bottom quarks, see Figs. [3] and @J 
The predictions at Q 2 = 50 GeV 2 and x = 10 -3 have been obtained with the 4-flavor PDF set of ETTl . 
Again, in the case of the MS mass the apparent convergence is much improved and, in particular, 
the relative size of NNLO a pp rox corrections is a few per cent only over the whole mass range 
considered for m^m^). 

In Figs. \5\ and [6] we investigate charm quark production in neutrino-nucleon DIS assuming 
an isoscalar target. We plot the CC charm structure function F^ for a nucleon at LO and NLO 
for Q 2 = lOGeV 2 , x = 10 _1 which corresponds to the typical kinematics of fixed-target neutrino- 
nucleon experiments. In comparison, the impact of higher order perturbative corrections is less 
than in the NC case discussed before. 

We compare the central predictions using a pole mass (Fig. [5] left) with those that employ 
the running mass (Fig. \5\ right) and the observed differences are rather marginal. This can easily 
be understood because the conversion to the MS mass only involves (e.g. for a W + -boson on 
isoscalar nucleon target) the derivatives of the up, down and the strange quark PDFs, which are 
all numerically rather small. Again, we also vary the renormalization and the factorization scale 
H r and fj. , independently by a factor of two around the central value \x 2 = n 2 = Q 2 + m 2 . The 
comparison is shown in Fig.[6]and the are no significant changes. In summary, we observe that the 
impact of the scheme change from a pole to a running mass is much more pronounced in the case 
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of NC DIS than in the case of CC DIS. 



We have demonstrated a clear improvement of the theoretical predictions for DIS heavy-quark 
structure functions by using the running mass. An immediate application of our results consists 
therefore in the direct determination of heavy-quark MS masses from the available DIS data. For 
reference in the following, we list the values for the charm and bottom masses in the MS scheme 
as quoted in the 2010 edition of the PDG [fTTTl . 

m c (m c ) = 1.27!°-°jGeV, (26) 



-0.09 

m b (m b ) = 4.19 ^GeV. (27) 



To start with, we can use the manifest dependence of the structure functions Fk on the heavy- 
quark mass to estimate the prospects of this approach. The relative uncertainty of such a mass 
determination is related to the corresponding uncertainty on the measurements of Fk as follows. 
Neglecting non-linear terms, a fit to the central prediction e.g. for F p in NC DIS (see Fig.Q]) yields, 

— ^ - 0.75-/. (28) 
m c F p 2 

This implies that a measurement of the proton structure function F p with an accuracy of 10% 
translates into a 0.75 x 10% = 7.5% uncertainty of the charm-quark mass. Thus, given the accuracy 
of current collider data from HERA (especially from HERA-II) an error on m c (m c ) of <9(few)% 
seems to be the ultimate precision one can aim at in this approach. 

For a quantitative comparison we conduct a phenomenological study similar to H19U29II , i.e. we 
perform a global fit of fixed-target (CCFR flU, NuTeV ED) and collider data 02113311 in the FFNS 
(with nf = 3) as a variant of ABKM ETTl . In the analysis we have taken the same 25 parameters 
as in [|2Ti which include also the strong coupling a s and the masses m c and nit, besides the usual 
PDF parameters. Interestingly, our fit does not return any sensitivity to the value of m^. Therefore 
we have constrained the bottom mass m^(m^) to its PDG value, i.e. Eq. (|27T) . For the running mass 
of the charm-quark, however, our analysis displays very good sensitivity and yields (depending on 
the order of perturbation theory) the following values 

m c (m c ) = 1.26 +0.09(exp) + 0.11 (th)GeV atNLO, (29) 

m c (m c ) = 1.01 ±0.09(exp) ±0.03(th)GeV at NNLO appr ox , (30) 

where the renormalization scale has been chosen n r = m c . We consider our mass determination 
at NNLOapprox accuracy as our best estimate. Eq. (|30l ) is the central result of this study and our 
determination is consistent with the world average at the level of +1.5cr for the quoted uncertain- 
ties. Our NNLO a pp r ox predictions are, of course, reliable only in a restricted kinematical range. 
However, given that they are generally rather small, we consider the agreement between the deter- 
minations at NLO and NNLO ap prox also a very good indication on the stability of the perturbative 
description. 

In Eqs. (|29l and (|3Q|) the experimental and theoretical uncertainties on m c (m c ) have been quoted 
separately. The former one is computed from the propagation of the statistical and systematic errors 
in the data, taking into account error correlations whenever available. The theoretical uncertainty 
is estimated from the sensitivity due to variations of the renormalization and factorization scales 



8 



ix r and fij. as follows. All current global PDF determinations assume = = Q in fits to DIS 
data (see e.g., the discussion in [26]). This is the appropriate scale choice for massless structure 
functions and, generally, for large values of Q when mass effects are negligible. In order to retain 
sensitivity to mass effects, especially in the region of low Q, we therefore determine the variation 
of F2 for the scale choice [i 2 = jj. 2 = Q 2 + Kin 2 , in the range k e [0, 8]. In this way, the variation 
of F2 at NLO (NNLO a pp r ox) results in the quoted uncertainty Am c (m c ) = +0.11 GeV (Am c (m c ) = 
+0.03 GeV). For consistency, we have also checked, that the statistical quality of our fit is not 
deteriorated, if we use these different scale choices, i.e. the obtained value for^ 2 changes by a few 
units only. 

In comparison to Refs. H211I34H the shapes of the PDFs and the value for the strong coupling do 
not change much in the present variant of the fit. We observe consistency within the +lcr statistical 
error. For illustration we display the + lcr band of absolute uncertainties for the (non-strange) light- 
quark (Fig. [7J left) and the gluon PDFs (Fig. [7J right) at the starting scale fi = 3 GeV compared to 
ABKM [21 J and we observe good agreement. The PDFs in Fig. [7J result from a fit, where we 
have additionally constrained the charm-quark mass by the PDG value of Eq. (|26l ). This results in 
m c (m c ) = 1.18 + 0.06(exp) + 0.03(th) at NNLO a pp r0 x in very good consistency with Eqs. ( 1261 and 
(1301) . For bottom the value m^imi,) = 4.19 + 0.12 with a symmetric error has been used. Since our 
analysis involves fixed-target data from CCFR/NuTeV, we have also paid particular attention to 
the strange-quark PDF and a potentially asymmetric strange sea. However, we find no indication. 
The total integrated asymmetry with the +lcr statistical uncertainty at the scale fi 2 = 20 GeV 2 is 
obtained as 



which is consistent with previous results [|29l . 

In Fig. [8] we confront data from NC heavy-quark DIS with the resulting predictions using 
running masses for charm and bottom, m c (m c ) = 1.18 and m^(m^) = 4.19. The ZEUS data [|32l 
displayed in Fig. [8] has not been used in the fit. At the smallest values of x and Q in Fig. [8] 
the predictions rise monotonically with increasing orders of perturbative QCD, thus improving 
agreement with the data. In this region the value of is sensitive to the coefficient functions for 
small ft, where the threshold approximation NNLO ap p r0 x is valid and can be considered as a good 
approximation to the full (yet unknown) NNLO result for F P 2 . As can be seen in Fig. [8l with the 
chosen value for charm (m c (m c ) = 1.18) our predictions are still slightly below data for at small 
values of x and Q. At large values of x and Q the slope of F~ flattens in Q, particularly for higher 
values of x, and the agreement with data is still very good. Future comparisons to high precision 
NC heavy-quark DIS data from the Run II of HERA will be interesting. 

A couple of interesting remarks can be made at this point. First of all, the experimental input 
to the PDG determinations of m c and originates entirely from e + e"-collisions or 5-decays (see 
e.g. 05113611 for recent analyses of e + e"-annihilation data with QCD sum rules). While the use of 
short distance masses is by now fairly standard in cross section predictions for those processes, it 
has not been used much in phenomenology at hadron colliders, although it is well known that the 
pole mass is plagued by large intrinsic ambiguities. However, with the increasing experimental 
precision of hadron collider data, there is a clear need to provide perturbatively stable theory pre- 
dictions and to use well-motivated definitions of fundamental quantities like the mass parameter. 

The MS masses of Eqs. ([291 and (1301) provide the first theoretically consistent determinations of 
these fundamental parameters in heavy-quark DIS, a process governed by the exchange of space- 




(31) 
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like bosons. So the agreement with the PDG value IfTTTl within the quoted accuracy is very reassur- 
ing also with respect to the different underlying (space-like) kinematics. Previously, charm mass 
determinations in heavy-quark DIS have been performed by CHARM II (3711 and NOMAD (381. 
They have extracted a value of m c from di-muon events in neutrino-nucleon DIS with a rather large 
uncertainty and based on a LO QCD analysis only. Also from the CCFR/NuTeV data for neutrino- 
nucleon DIS a value for the charm mass has been extracted (391 , which however is effectively LO 
only, as far as the mass dependence of the cross section is concerned. None of these determinations 
enters the world average as quoted by the PDG. 

Another interesting issue concerns the heavy-quark PDFs. These are needed at high-energy 
hadron colliders for hard scattering processes at scales Q :» m c ,mb, where a FFNS with effectively 
nf = A ox nf = 5 light flavors is the appropriate description. The PDFs for charm- and bottom- 
quarks in 4- and 5-flavor schemes can be generated from the ones obtained in a 3-flavor FFNS as 
convolutions of the gluon and flavor-singlet distributions with massive operator matrix elements. 
Through the explicit mass dependence of the latter the uncertainty on heavy-quark PDFs is directly 
related to the accuracy of the quark mass parameter. Preliminary studies have shown that precision 
input for the values of charm and bottom masses can greatly improve the accuracy of charm- and 
bottom-quark PDFs. Future studies will be devoted to an in-depth investigation of the transitions 
from 3- to 4- to 5-flavor FFNS with the running mass scheme (40l . 

Finally, let us address the implications for LHC processes, which are quite clear. In global PDF 
fits, the predicted rate for W ± - and Z-boson production is very sensitive to the chosen pole mass 
value for charm in schemes with four or five active light flavors. At the = 7 TeV LHC, for 
instance, Ref. (4Tf. has reported shifts in the total cross sections for W ± - and Z-boson production 
of more than 2% resulting from small variations of 0.15 GeV in the pole mass value for m c around 
a central value (see Tab. [I]). For the VS 7 " = 14 TeV LHC, these uncertainties increase to more than 
3%. These findings are worrisome as they potentially invalidate precision predictions for these 
important high precision measurements of Standard Model processes. 

Fortunately, as we have demonstrated, these uncertainties can be almost entirely eliminated 
by adopting the running mass. By using e.g. the very precise world averages as constraints in 
global fits, the errors on m c and can be greatly reduced. Moreover, these errors can be directly 
propagated to the PDF uncertainties with no need for additional assumptions. In this manner, very 
accurate and precise predictions for W ± - and Z-boson production at the LHC are possible. 

4 Summary 

We have studied the production of heavy quarks in NC and CC DIS including higher order radiative 
corrections in QCD adopting the running MS mass and we have demonstrated the clear advantage 
of using this scheme. The resulting theory predictions display an improved apparent convergence 
through NNLOapprox as compared to the result in the pole mass scheme, especially in the NC case. 
Also the stability of the perturbative series with respect to scale variations is much greater. 

We have used our results to determine for the first time the MS mass for the charm-quark m c (m c ) 
based on first principles in QCD from a fit to DIS data for heavy-quark production. The obtained 
value is consistent with the world average as published by the PDG and it provides complementary 
information on this fundamental parameter from hadronic processes with space-like kinematics. 

We have shown that the use of MS masses in heavy-quark DIS can also improve predictions 
for hadron colliders by eliminating sizable uncertainties in PDFs. This underpins the need for 
global fits of PDFs to adopt the running mass scheme. Current global fits of PDFs employ the pole 
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mass scheme and small variations in the chosen value for the charm mass can easily amount to 
differences of a few per cent in the predicted cross sections of W ± - and Z-bosons at LHC energies. 
With a physically well motivated short distance mass these uncertainties can be eliminated to a 
large extent. Progress in this direction will be reported elsewhere 11401 . 

The numerical code for the computation of structure functions for heavy-quark production in 
deep-inelastic scattering with a running mass is publicly available for download from 11421 or from 
the authors upon request. 
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Figure 1: The mass dependence of the NC charm structure function F p for a proton with Q 2 = 10 GeV 2 , 
x = 10~ 3 and \x r = \Xj- V<2 2 + 4ra 2 using the PDFs of (2TJ. The charm-quark mass is taken in the on-shell 
scheme (left) and in the MS scheme (right) at LO (blue), NLO (green) and NNLO a p pr ox (red). 
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Figure 2: Same as in FigQ] The band denotes the independent variation of the scales n r ,^i^ - k V<2 2 + 4m 2 
in the range k e [1 /2, 2]. 
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Figure 3: The mass dependence of the NC bottom structure function F P 2 for a proton with Q 2 = 50 GeV 2 , 
x = 10~ 3 and \i r = l^f- ^Q 2 + 4m 2 . The bottom-quark mass is taken in the on-shell scheme (left) and in the 
MS scheme (right) at LO (blue), NLO (green) and NNLO approx (red). 
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in the range k e [1/2,2]. 



_ 1 1 1 1 1 1 1 1 1 

: p 2 charm (ABKM 3 NLO) 


i i i I 
LO 


1 1 


rQ 2 =10GeV 2 , x = 10" 1 


NLO 




L u 2 = Q 2 + m c 2 






: , , I , , , , I , 


, , , I 





1 1.5 2 

pole mass m (GeV) 



0.3 
0.25 

0.2 
0.15 

0.1 
0.05 




- I 1 1 1— 



~ I 1 1 r~ 



F 2 charm (ABKM 3 NLO) LO 
Q 2 = 10GeV 2 , x=10" 1 NLO 
H 2 = Q 2 + m 2 




_L 



_l I I L_ 



_L 



_l I I l_ 



1 1.5 2 

running mass m (m ) (GeV) 



Figure 5: The mass dependence of the CC charm structure function for a nucleon with Q 2 = 10 GeV 2 , 
x = 10 -1 and [x r = fij. = -\JQ 2 + m 2 using the PDFs of [21]. The charm-quark mass is taken in the on-shell 
scheme (left) and in the MS scheme (right) at LO (blue) and NLO (green). 
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Figure 6: Same as in Figj5] The band denotes the independent variation of the scales fj. r ,fi, = k^Q 2 + m 2 
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H=3 GeV 




Figure 7: The light-quark (left) and gluon (right) PDFs obtained in the global fit: The dotted (red) lines 
denote the ±lcr band of absolute uncertainties resulting from the fit of this paper and the solid (red) line 
indicates the central prediction with with running masses 1.18 GeV (charm) and 4.19 GeV (bottom). For 
comparison the shaded (grey) area represents the results of ABKM [21 ]. 
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Figure 8: The predictions of the global fit for the charm structure function F L 2 compared to data of Ref. ||32| 
on F L 2 from Run I of HERA. The running mass m c (m c ) = 1.18 GeV has been obtained in the variant of the 
fit with the PDG value Eq. (l26l) as an additional constraint. 
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